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Abstract. Two peculiar stars showing an apparent extremely broadened and 
strong NalD absorption have been discovered in surveys for cool white dwarfs 
by Oppenheimer et al. (2001) and Harris et al. (SDSS, 2003). We discuss 
the nature of these objects using PHOENIX atmosphere models for metal-poor 
brown dwarfs/very low mass stars, and new white dwarf LTE and NLTE models 
for hydrogen- and helium-dominated atmospheres with metals. These include 
complete molecular formation in chemical equilibrium and a model for the al- 
kali resonance line broadening based on the damping profiles of Allard et al. 
(2003), as well as new molecular line opacities for metal hydrides. First results 
of our calculations indicate good agreement with a hydrogen-dominated WD 
atmosphere with a Na abundance roughly consistent with a state of high ac- 
cretion. We analyse deviations of the abundances of Na, K, Mg and Ca from 
the cosmic pattern and comment on implications of these results for standard 
accretion scenarios. 
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Introduction 



WD2356— 209 is an apparent cool white dwarf id entified in the high prope r 
motion-survey for halo white dwarf candidates of lOppenheimer' et all (l2nnih . 
with a spectrum showing very unusual wide and d eep absorpti o n at 5000 - 
6000 A. Its exotic properties have been confirmed bv ISalim et alJ ( 20041 ) . who 



found a combination of blue B — V and extremely red V — I colours putting 
this star far outside of observed and theoretical sequences of both hydrogen- 
and helium-atmosphere white dwarfs. Meanwhile, another star with a similar 
spectroscopic and photometric appe arance had been identified in Data Releas e 
2 of the Sloan Digital Sky Survey (iHarris et alJ '200.'^': ' Kleinman et al.ll2004 ): 
SDSS J 13 30 01.13+64 35 23.8 (hereafter SDSS J1330+6435), which is also show- 
ing high proper motion of O'.'193/yr. The prominent flux depression in V in these 
stars has been suggested to be due to a strong and heavily broadened Na I D dou- 
blet. Similar features have been observed in the coolest brown dwarfs of spectral 
type T and identified with the resonance lines of the alkali metals Na and 
broadened by van d er Waals interaction with both H2 and He ( Tsuii et al]ll999l : 



iBurrows et ai]l200nh . These observations made clear that at line widths of more 
than some 10 A the standard impact theory of line broadening, leading to a clas- 
sical Lorentz profile of the far wings, could no longer adequately reproduce the 
observed line shapes. More detailed quantum-mechanical calculations of the in- 
teraction between alkali atoms and H2 and He perturbers ( Burrows Sz VolobuvevI 



I2OO3I : lAllard et al.ll200^ have shown that the Lorentzian line profile underesti- 
mates the absorption cross section at distances beyond 2-3 FWHM from the 
line centre, out to ~ 2000 A, while at very large detuning it significantly over- 
estimates it. In this paper we are exploring theoretical spectra for a number of 
different stellar models, using our latest calculations of alkali line profiles, and 
evaluate their capability to explain the observed features of WD2356— 209 and 
SDSSJ1330+6435. 



2. Model Atmospheres 

We hav e calculated mode ls using the general stellar atmosphere code PHOENIX 



V. 13.6 ((Hauschildt Sz Baron..l999. 1. PHOENIX has been continuously developed to 



model the atmospheres of the coolest stars, brown dwarfs and gas giant planets, 
incorporating an equation of state (EOS) that treats over 600 molecular species 
in chemical equilibrium (CE) and includes the formation of liq uid and solid 
condensates ( AUard et al.ll200lh . We are using the unified theory of lAllard et al.l 



(|2003f ) to predict the complete spectrum of the resonance lines of the alkali metals 
Li, Na, K, Rb and Cs from the centre to the extreme wings, however with 
normal abundance patterns only Na and K are found to have widths beyond 
some tens of A, and only for these perturbations due to both neutral He and H2 
in two interaction geometries are available at this time. These models have been 
successfully applied to the highly transparent and dense atmospheres of cool T 
dwarfs, the only astrophysical objects in which such lines have been observe d 
prior to the two stars that are the topic of this paper (lAllard et al.ll2003l . l2004^ . 
The colours of WD2356-209 and the blue spectrum of SDSS J1330+6435 suggest 
that the A 4227 A Cal line is also an important absorber, but for this earth 
alkali metal no detailed calculations of the line spectrum are available yet. We 
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have tried to simulate the broadening of this hne by using the Nal data, but 
multiplying the damping constant with GO*^'^ to account for its non-hydrogenic 
configuration. 

In other respects though, T dwarf spectra disagree very strongly with the 
observed signatures of WD2356— 209 and SDSS J1330+6435: they show equally 
prominent absorption from the K I AA 7667,7701 A resonance doublet as from Na, 
and they have extremely red optical-to-near IR spectra, with the optical flux al- 
most completely quenched bluewards of the NalD doublet. Objects this cool 
could therefore not account for the quite blue B — V colours of WD2356— 209. 
For any brown dwarf or very low mass-star of higher temperature (Teg ^ 1400 K) 
however the strength of alkali line absorption rapidly decreases, as the atmo- 
spheres become more opaque due to strong absorption from molecular bands 
and condensates. Keeping a strong signature of the resonance lines in atmo- 
spheres of higher temperature would require a metallicity significantly below 
the solar value. The thusly reduced amount of molecules and dust and would 
increase the transparency of the atmosphere, allowing the atomic lines to form 
again under conditions of high density and pressure. To explore this possibility 
we have therefore calculated a grid of dwarf atmosphere models at Tgfj up to 
3000 K and [Fe/H] down to —5. For the main alternative explanation, these 
stars indeed being white dwarfs with a high abundance of photospheric metals, 
we computed models with log 5 of 7.0-8.5 and hydrogen- or helium-dominated 
atmospheres . For a ll models the partial pressures of all molecules in the EOS of 
lAUard et al.1 (l200lh have been calculated in full CE, although in the hydrogen- 
rich WD models chemistry is dominated by only a few hydrogen compounds, 
and in helium-rich ones of course not many molecules form at all. The models do 
not include, on the other hand, the effects of a non-ideal EOS that may become 
important in white dwarfs this cool and dense (cf. also Kowalski, Mazevet & 
Saumon, this vol.!). 



3. Comparison vi^ith Models 
3.1. The Nature of the Dwarf (s) 

Figure^ shows a representative spectrum for an extremely metal-poor model 
in the L-dwarf temperature range, as might be conceived for a Population III 
star just above the hydrogen-burning limit. It is immediately apparent that, 
without the dust opacity dominating normal L-dwarf spectra, the alkali lines 
can be strong enough to qualitatively explain the observed absorption features. 
However the line profiles fail to reproduce the observed shapes in detail, and 
show still too red optical— near-IR colours. We conclude that a low mass-star or 
brown dwarf- type of model can be excluded for these two stars. This makes the 
most likely classification of these objects white dwarfs with traces of metal in an 
either helium- (DZ) or hydrogen-dominated atmosphere (we have denoted the 
latter as DAZ, in analogy to other hydrogen-rich white dwarfs showing metal 
absorption lines besides the Balmer series of hydrogen, although due to the 
extremely low temperatures in this case, Balmer lines would probably not be 
observable even at the highest resolutions). 



4 



Homeier et al. 



H-models log g . 8.0, [M/HJ.-3.5, [Na,K/H]=-1 .5 




Figure 1. a: (left) Comparison of a [Fe/H] = — 4 VLM-star model to the 
spectra of WD2356-209 and SDSS J1330+6435. The SDSS spectrum has 
been averaged with a 4-pixel boxcar filter to i? w 1000. b: (right) Fit of 
hydrogen-rich WD-models to the spectrum of SDSS J1330-I-6435. Principal 
absorption features are indicated 



3.2. White Dwarf Models 

Hydrogen- dominated atmospheres A series of fits to WD2356— 209 to hydrogen- 
LTE atmospheres is presented in Fig. ^p. The strong dependence of the alkali 
lines on temperature, becoming quickly weaker with increasing Tefj, is evident. 
Test models for solar abundance patterns predicted the formation of stronger 
molecular bands, chiefly due to MgH and CaH, than are present in the obser- 
vations. Although the spectrum of WD2356— 209 shows an apparent additional 
opacity source in the blue wing of the Nal doublet, its position does not quite 
agree with the modelled band, while SDSS J1330-I-6435 only shows a marginally 
detectable feature. In addition, CaH bands should be visible around 7000 A if 
the elemental Ca abundance were in the same ratio to Na as in the Sun. The 
models shown here therefore have Na and K abundances set to — l.Sdex relative 
to solar, and all other metals to —3.5. As can be seen, even at T^s < 4000 K 
quite high Na abundances are required to reproduce the strength of the D lines. 
In contrast the KI doublet, which is marginally if at all detectable in the obser- 
vations, is predicted to be much stronger at that abundance. [K/H] therefore 
seems to have a lower abundance more similar to the other metals. 



Helium-rich models In addition to the degeneracy of solutions with respect 
to Teff and metal abundances, the main constituent in the atmosphere is also 
unknown. The effect of a varying H/He mixture is illustrated in Fig. [2K- We 
find a slight improvement of the line fit with moderate amounts of helium, e. g. 
50-90%. Infrared flux is also more suppressed, due to the stronger H2-He 
collision induced absorption (CIA), with increasing helium abundance even up 
to H/He = 10~^^. These models thus show rather blue synthetic IR colours of 
V — H^ 2.0, as compared 3.0-3.5 for pure hydrogen models and an observed 
V-H = 2.75 for WD2356-209 (Bergeron et al, this vol.). The best fits seem to 
be obtained with H/He of order unity. We did not run any completely hydrogen- 
free models, which should eventually remove the CIA beyond 1 fim. Bergeron 
et al. report a good fit with such models, although they only use standard 
Lorentzian broadening for the Nal doublet. 
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Ts„ . 3500K log g = 8.0 H-models Ten = 3800K log g . 8.0, [M;H]=-3.5, [Na.K/H]=-1 .5 




Figure 2. a: (left) Comparison of models with different H/He number ratios 
for SDSS J1330+6435. b: (right) NLTE effects on H-atmosphere models. High 
departure ratios in the upper atmosphere may be responsible for narrower line 
cores (long-dashed plot), however coUisional deexcitation due to H2 would 
most likely restore LTE (cf. text) 



NLTE Effects The line profiles particularly in our pure-hydrogen LTE models 
appear generally too deep in the near wings out to a few 100 A, indicating that 
either the far wing is too weak compared to the core and near wing in our line 
profiles, or for some reason Nal absorption is overestimated in the upper atmo- 
sphere, where the inner parts of the lines form. NLTE effects could contribute 
to the latter. Departures from LTE are typically observed in white dwarfs of 
Tgfj ^ 20 000 K, and might seem unlikely in dwarfs this cool and dim. How- 
ever, the highly neutral atmosphere of ultracool dwarfs leaves almost no free 
electrons as the main source of collisional deexcitation and could thus favour 
NL TE effects. Alkali line absorption under such circumstances has been studied 
by iBarman etaP (|2002l ) for the transit planet HD 209458b. Due to the ex- 
treme paucity of free electrons the dominant species, molecular hydrogen, will 
most likely become the driving force to restore LTE. This effect can not be pre- 
cisely modelle d, since no cross-sections for H2-Na collisions are available yet. 
iBarman et al.l have discussed two limiting cases, including only e~ collisions, 
and including H2 collisions with the same cross-sections as for e~. Theory sug- 
gests a much smaller H2 cross-section, scaling with the inverse mass ratio. We 
have therefore also calculated NLTE models using 4000 and 40 000 times smaller 
cross-sections for H2 . The test models in Fig. with e~ only show very signif- 
icant deviations from the LTE line shape, but even with the lowest estimate for 
H2 collisional rates LTE is almost completely restored. NLTE effects are thus 
probably not the cause of the discrepancies between models and observation. 



Limitations of the models Since LTE does seem to hold within the line- forming 
regions, actual depletion of elemental Na by gravitative diffusion could be an- 
other explanation for the missing opacity in the deep line cores. This is not 
an uncommon effect in white dwarf atmospheres, but a detailed calculation of 
it is beyond the scope of our current models. We also have to caution, how- 
ever, that the calculations for the line spectra we used were designed to be used 
with perturber densities of up to 10^^ cm~^, whereas at the higher densities 
of white dwarf atm ospheres the single-perturber approximation breaks down 
( Allard et all 1200.^ . heading to an underestimation of the far wing part of the 
line. 
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4. Conclusions 

With the low resolution and S/N of both observations and the uncertainties 
of our current models it is not possible to draw strong conclusions about the 
atmospheric conditions in these white dwarfs. E. g. the presence of hydride 
bands cannot be clearly established, which would allow a better determination 
of elemental abundance ratios, and confirm whether hydrogen is an important 
constituent of the atmosphere. The quite high Sodium abundance required in the 
DAZ scenario, which might be beyond what interstellar accretion could provide, 
is an obvious argument against such models. The apparent overabundance of Na 
relative to Ca, and perhaps also Mg and K, also disagrees with selective dust- 
accretion models, which predict rather lower abundances of the volatile Alkalis 
((Zuckerman et al. 2003). However, DZ models would face the same irregular 
abundance patterns, except that hydride bands are not expected. Also the 
problem of avoiding the accretion of hydrogen from the ISM is still unresolved. 
With the current knowledge, a model with He/H ratio of order unity might best 
reconcile both observed features and the problems with the extremely short 
diffusion timescales of pure hydrogen atmospheres. 

Observations that could resolve molecular bands, and possibly lines from 
higher excitation and ionisation levels for better constraints on temperature, 
would be of great help in deciding between these different interpretations. New 
theoretical work on the behaviour of atomic lines at very high perturber densities 
should also both benefit from and aid in the study of these unusual white dwarfs. 
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